FRACTAL CONTINUATION 



MICHAEL F. BARNSLEY AND ANDREW VINCE 

Abstract. A fractal function is a function whose graph is the attractor of an 
iterated function system. This paper generalizes analytic continuation of an 
analytic function to continuation of a fractal function. 



1. Introduction 

Analytic continuation is a central concept of mathematics. Riemannian geometry 
emerged from the continuation of real analytic functions. This paper generalizes 
analytic continuation of an analytic function to continuation of a fractal function. 
By fractal function, we mean basically a function whose graph is the attractor of an 
iterated function system. We demonstrate how analytic continuation of a function, 
defined locally by means of a Taylor series expansion, generalises to continuation 
of a, not necessarily analytic, fractal function. 

Fractal functions have a long history, see [TS] and [HI Chapter 5]. They were 
introduced, in the form considered here, in pQ. They include many well-known types 
of non-differentiable functions, including Takagi curves, Kiesswetter curves, Koch 
curves, space-filling curves, and nowhere differentiable functions of Weierstrass. A 
fractal function is a continuous function that maps a compact interval / cR into 
a complete metric space, usually K or R 2 , and may interpolate specified data and 
have specified non-integer Minkowski dimension. Fractal functions are the basis 
of a constructive approximation theory for non-differentiable functions. They have 
been developed both in theory and applications by many authors, see for example 
HI H3 HH 02 Q31 QI] and references therein. 

Let N be an integer and 1 = {1, 2, N}. Let M > 2 be an integer and X C R M 
complete with respect to a metric dx that induces, on X, the same topology as the 
Euclidean metric. Let W be an iterated function system (IFS) of the form 

(1) W = {X; W „,neI}, 

We say that W is an analytic IFS if w n is a homeomorphism from X onto X for 
all n£l, and w n and its inverse if" 1 are analytic. By w n analytic, we mean that 

w n (x) = (w nl (x),w n2 (x),...,w nM (x)), 

where each real- valued function w nm (x) — w nm (xi,X2, im) is infinitely differen- 
tiable in Xi with xj fixed for all j ^ i, with a convergent multivariable Taylor series 
expansion convergent in a neigbourhood of each point (xi,x%, ...,Xm) G X. 

To introduce the main ideas, define a fractal function as a continuous function 
/:/—>• M M_1 , where / C M is a compact interval whose graph G(f) is the attractor 
of a IFS for the form in eqaution p~| . A slightly more restrictive definition will be 
given in Section [3] If W is an analytic IFS, then / is called an analytic fractal 
function. 
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The adjective "fractal" is used to emphasize that G(f) may have noninteger 
Hausdorff and Minkowski dimensions. But / may be many times differentiable or 
/ may even be a real analytic function. Indeed, we prove that all real analytic 



functions are, locally, analytic fractal functions; see Theorem 4.2 An alternative 



name for a fractal function / could be a "self-similar function" because G(f) is a 
union of transformed "copies" of itself, specifically 

N 

(2) G(f) = |J w n (G(f)). 

n=l 

The goal of this paper is to introduce a new method of analytic continuation, 
a method that applies to fractal functions as well as analytic functions. We call 
this method fractal continuation. When fractal continuation is applied to a locally 
defined real analytic function, it yields the standard analytic continuation. When 
fractal continuation is applied to a fractal function /, a set of continuations is 
obtained. We prove that, in the generic situation with M = N = 2, this set of 
continuations depends only on the function / and is independent of the particular 
IFS W that was used to produce /. The proof relies on the detailed geometrical 
structure of analytic fractal functions and on the Weierstass preparation theorem. 



Figure 1 . This paper concerns analytic continuation, not only of 
analytic functions, but also of non-differentiablc functions such as 
the one whose graph is illustrated here. 

The spirit of this paper is summarized in Figure [T] Basic terminology and back- 
ground results related to iterated function systems appear in Section [2] In Section 
[3] we establish the existence of fractal functions whose graphs are the attractors of a 
general class of IFS, which we call an interpolation IFS. An analytic fractal function 
is a fractal function whose graph is the attractor of an analytic interpolation IFS. 
This includes the popular case of affine fractal interpolation functions pQ . An ana- 
lytic function is a special case of an analytic fractal function, as proved in Section [4] 
Fractal continuation, the main topic of this paper, is introduced in Section [5] The 
fractal continuation of an analytic function is the usual analytic continuation. In 
general, however, a fractal function defined on a compact domain, has infinitely 
many continuations, this set of continuations having a fascinating geometric struc- 
ture as demonstrated by the examples that are also contained in Section [5] The 
graph of a given fractal function can be the attractor of many distinct analytic 
IFSs. We conjecture that the set of fractal continuations of a function whose graph 
is the attractor of an analytic interpolation IFS is independent of the particular 
IFS. Some cases of this uniqueness result are proved in Section [6j 
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2. Iterated Function Systems 
An iterated function system (IFS) 

W = {X;w n ,nel} 

consists of a complete metric space X C R M with metric dx, and N continuous 
functions w n : X — > X. The IFS W is called contractive if each function w in W 
is a contraction, i.e., if there is a constant s G [0, 1) such that 

d(w(x),w(y)) < sd(x,y) 

for all x, y G X. The IFS W is called an invertible IFS if each function in W is 
a homeomorphism of X onto X. The definition of analytic IFS is as given in the 
introduction. The IFS W is called an affine IFS if X = R M and w n : R M -> R u 
is an invertible affine map for all n £ I. Clearly an affine IFS is analytic, and an 
analytic IFS is invertible. 

The set of nonempty compact subsets of X is denoted H = H(X). It is well- 
known that H is complete with respect to the Hausdorff metric h, defined for all 
S, T G H, by 

h(S,T) — max < max mm dx(s.t), max mindx(s,t) 
I ses ter ter ses 

Define W : H -> H by 

(3) W{K) = |J w n (K) 

nel 

for all K G H. Let W° : H ->• H be the identity map, and let W k : i-> H be the 
fc-fold composition of W with itself, for all integers fe > 0. 

Definition 2.1. A set A G H is said to be an attractor of W if W(A) = A, and 

(4) lim W fe (A') = A 

for all AT G H, where the convergence is with respect to the Hausdorff metric. 

A basic result in the subject is the following [7]. 
Theorem 2.2. If IFS W is contractive, then W has a unique attractor. 

The remainder of this section provides the definition of a certain type of IFS 
whose attractor is the graph of a function. We call this type of IFS an interpolation 
IFS. We mainly follow the notation and ideas from [2 [21 E] . 

Let M > 2 and N an integer and 1= {1,2, N}. For a sequence xq < 
xi < ■■■ < x j\[ of real numbers, let L n : K — > K be the affine function and 
F n : K x M M_1 — > _R M_1 , n = 1,2,..., A be a continuous function satisfying 
the following properties: 

(a) L n (x ) = x n -i and L n (x N ) = x n . 

(b) There are points yo and yw in IR M_1 such that Fi(xo,yo) = Ho an d 
F N (x N ,y N ) = y N . 

(c) F„ + i(a;o,yo) = F n (x N , y N ) for n = 1, 2, . . . , N - 1. 
Let W be the IFS 

(5) W = {M M ; w n , n G I}, 
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where 

(6) w n (x, y) = (L n (x),F n (x, y)). 
Keeping condition (c) in mind, if we define, for each n£l, 

y n := F n+1 (x Q ,y ) = F n (x N ,y N ), 

then note that 

(7) w n {x ,yo) = {x n -\,Vn-\) and w n (x N , y N ) = (x n , y n ). 

Definition 2.3. An interpolation IFS is an IFS of the form given by ( |5~| and 
( |6~| above that satisfies (1) w n is a homeomorphism onto its image for all n G I, 
and (2) there is an s G [0, 1) and an A4 G [0, oo) such that 

(8) \F n (x,y)-F n (x',y')\<M\x-x'\+s\y-y'\ 

for all x,x' G M.,y,y' G ]R^ /_1 and for all n E I. The term "interpolation" is 



justified by statement (2) in Theorem 3.2 in the next section. 



3. Fractal Functions and Interpolation 



Properties of an interpolation IFS are discussed in this section. Theorem 3.2 
the main result. 

Lemma 3.1. If W is an interpolation IFS, then W is contractive with respect to 
a metric inducing the same topology as the Euclidean metric on M. M . 

Proof. Let L n — a n x + b n and let d be the metric on R M defined by 

d((x,y),(x',y')) = e\x-x'\ + \y-y'\ 

where e G (A4/(l— a), oo) and a = max{a„ : n € I}. The metric d is a version of the 
"taxi-cab" metric and is well-known to induce the usual topology on R A1 . Moreover, 
w n is a contraction with respect to the metric d: for (x, y), (x' , y') G K x M. x , 

d(w n {x,y),w n (x',y')) = e\L(x) - L(x')\ + \F n (x,y) - F n {x' ,y')\ 

< ea\x - x'\ + M\x - x'\ + s\y - y'\ 
= (ea + M) \x - x'\ + s\y - y'\ 

— ce\x — x'\ + s\y — y'\ 

< max{c, s} d((x,y), (x',y')) 

where c = a + A4/e is monotone strictly decreasing function of e for e > 0, so c is 
strictly less than a + M/(M/(1 - a)) = 1 for e G {M/(l - a),oo). □ 



Theorem |3 . 2 1 generalizes results such as [SI p. 186, Theorem 5.4]. Let I = [x<j, a; at], 
C(I) = {/:/—> M M_1 : / is continuous} and C (I) := {g G C{I) : g(x ) = 

yo,g{x N ) = y N }. 

Theorem 3.2. IfW is an interpolation IFS, then 

(1) The IFS W has a unique attractor G :— G(f) that is the graph of a con- 
tinuous function f : I — > 

(2) The function f interpolates the data points (xo, yo), (xx, y±), . . . , (x^, J/iv)> 
i.e., f(x n ) = y n for all n. 
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(3) If W : Coil) -> Co(I) ™ defined by (Wg)(x) = F n {L~ l [x) , g{L~ x {x))) for 
x £ [x n _i, x n ], for n€l, and /or all g £ Cq{I), then W has a unique fixed 
point f and 

f= km W k (f a ) 

n— >oo 

for any f £ C Q (I). 

Proof. It is readily checked that the mapping W of equation p~| takes Cq(I) into 



Cq(/) and also that the mapping W of statement (3) in Theorem 3.2 takes Cq{I) into 



C (J). Moreover, if G(/ ) is the graph of the function / G C (I), then W(G(/ )) 



W-^/o). This implies, by property (4 ) of the attractor, that the function / in 
statement (1), assuming that it exists, is the same as the function / of statement 
(3), assuming that it exists. Statement (3) is proved first. 

(3): That the map W is a contraction on Co(I) with respect to the sup norm 
can be seen as follows. For all g £ Co (J), 

\(W 9l )(x) - (Wg 2 )(x)\ < m a x\F n (L- 1 (x),g 1 (L- 1 (x))) - F n (L~ 1 (x),g2(L~ 1 (x)))\ 

< maxs \giiL-\x)) - g 2 (L~ 1 (x))\ 

< s \gi(x) - g 2 (x)\ , 

where < s < 1 is the constant in condition ( |8""| ). Statement (3) now follows from 
the Banach contraction mapping theorem. 



(1) : According to Lemma 3.1 the IFS W is contractive. By Theorem 2.2 W has 
a unique attractor G. Let Go denote the graph of some function f in Co(I). Using 
statement (3) there is a function / £ Cq(I) such that / = lim/^oo W k (fo). By 
what was stated in the first paragraph of this proof and by the property p~| in the 
definition of attractor, we have G = lim^oo W fc (Go) = G(limfe_ J . 00 W fe (/o)) = <?(/)■ 

(2) : The attractor G must include the points (xo,yo) and (xn, J/w ) because 
they are fixed points of W\ and w^. Hence, by Equation ( |7 [ ) G must contain 
(x n ,y n ) = w n (x N ,y N ) for all n. □ 



Remark 3.3. Theorem 3.2 remains true if F„ : X — > R. and w„ : X — > X, for 



all n £ I, where X C M. M is a complete subspace of K M , and X contains the line 
segment [a;o,XAr] (treated as a subset of E M ). This, for example, is the situation in 
Theorem IL2I in the next section. 



Definition 3.4. A function / whose graph is the attractor of an interpolation IFS 
will be called a fractal function. A function / whose graph is the attractor of an 
analytic interpolation IFS will be called an analytic fractal function. Note that, 
although a fractal function usually has the properties associated with a fractal set, 
there are smooth cases. See the examples that follow. 

Example 3.5. (parabola) The attractor of the affine IFS 

W={R 2 ; Wl (x,y) = {x/2,y/A),w 2 (x,y) = ((x + l)/2, (2x + y + l)/4)} 

is the graph G of / : [0, 1] — > M, f(x) = x 2 . For each k £ N the set of points 
W fc ({(0,0),(l,l)}) is contained in G and the sequence {W k ({{0, 0), (1, 1)})}£L 
converges to G in the Hausdorff metric. Also, if /o : [0, 1] is a piecewise affine 
function that interpolates the data {(0,0), (1, 1)} then f k := W k (f ) is a piecewise 
affine function that interpolates the data W fe ({(0, 0), (1, 1)}) and {fk}kL Q converges 
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to / in (C([0, 1], doo). The change of coordinates (x, y) — > (y, x) yields an affine IFS 
whose attractor is the graph of / : [0, 1] — > R, f(x) = y/x. 

Example 3.6. (arc of infinite length) Let d\ + c?2 > 1, di,d 2 & (0,1)- The 
attractor of the affine IFS W = {M 2 ; w 2 }, where 
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is the graph of / : [0,2] R which interpolates the data (0,0), (a, c), (1,0) and 
has Minkowski dimension D > 1 where (see e.g. [9j p. 204, Theorem 5.32]) 

(a) D - 1 d 1 + {l-a) D - 1 d 2 = l. 

Example 3.7. (once differentiable function) The attractor of the affine IFS 
W = {R 2 ; wi, w 2 }, where 
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is the graph of a once differentiable function / : [0, 2] — > E which interpolates the 
data (0,0), (2/3,1), (2,0). The derivative is not continuous. The technique for 
proving that the attractor of this IFS is differentiable is described in [5]. 

Example 3.8. (once continuously differentiable function) The attractor of 
the affine IFS W = {R 2 ;w 1 ,w 2 }, where 
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is the graph of a continuous function / : [0, 1] — ► M. that possesses a continuous first 
derivative but is not twice differentiable, see [H[S]. 

Example 3.9. (nowhere differentiable function of Weierstrass) The attrac- 
tor of the analytic IFS 

W = {R 2 ;w 1 (x,y) = (x/2, £y + sinTra;), w 2 (x, y) = {{x + l)/2, £y - sinTrx)} 

is the graph of / : [0, 1] — >• M well-defined, for |^| < 1, by 

00 

(9) f(x) =^2$ k sin2' £+1 7ra; for all x € [0, 1], 

fc=0 

This function f(x) is not differentiable at any x £ [0, 1], for any £ 6 [0.5, 1). See 
for example [51 Ch. 5]. 

4. Analytic Functions are Fractal Functions 

Given any analytic function / : I R, we can find an analytic interpolation 
IFS, defined on a neighbourhood G of the graph G(f) of /, whose attractor is G(f). 
This is proved in two steps. First we show that, if f'(x) is non zero and does not 
vary too much over /, then a suitable IFS can be obtained explicitly. Then, with 
the aid of an affine change of coordinates, we construct an IFS for the general case. 
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Lemma 4.1. Let f : [0, 1] — > K fee analytic and strictly monotone on [0,1], 
urei/i bounded derivative such that both (i) max xg r 01 i |/'(x/2)//'(x)| < 2 and (u,) 
max^grQ^i |/'((x + l)/2) /f'(x)\ < 2. T/ien i/iere is a neighbourhood Get 2 of 
G(f), complete with respect to the Euclidean metric on M 2 , such that 

W ={G; Wl (x, y) = (x/2, f(f- 1 {y)/2)),w 2 (x, y) = (x/2 + 1/2, /(/- 1 (y)/2 + 1/2))} 

is an analytic interpolation IFS whose attractor is G(f). 

Proof. First note that G — G(/) is compact and nonempty. Also G is invariant 
under W because 

W(G)=/i(G)U/ 2 (G) 

= {(x/2J(f- 1 (f(x))/2):x£ [0,1]}U 

{(x/2 + 1/2, f{f-\f(x))/2 + 1/2) : x € [0, 1]} 

= {(x/2, /(x/2)) : x € [0, 1]} U {(x/2 + 1/2, /(x/2 + 1/2) : x E [0, 1]} 

= G. 



Second, we show that property (|8j) holds on a closed neighbourhood G of G. To 
do this we (a) show that it holds on G and then (b) invoke analytic continuation 
to get the result on a neighbourhood of G. Statement (a) follows from the chain 
rule for differentiation: for all y such that (x, y) £ G we have, 



dy 



f{f-\y)m 



f {f-\y)m 



2/'(/-Hy)) 



f'(x/2) 



2f'(x) 



(since y = /(x) on G) 



< 1. 



To prove (b) we observe that both v\(y) = f{f~ 1 (y)/2) and vi = f(f~ 1 (y)/2 + 
1/2) are contractions for all y such that (x,y) e G and so, since they are both 
analytic, they are contractions for all y in a neighbourhood N of /([0, 1]). Clearly 
x/2 and x/2 +1/2 are contractions for all x G K. Finally, it remains to show that 
G cJ\f can be chosen so that w n (G) cG (n = 1,2). Let G be the union of all 
closed balls B(x, y) of radius e > 0, centered on (x, y) e G, where e is chosen small 
enough that G cAf. 

The inverse of wj -1 : (G) -» G is well-defined by w^ix^) = (2x, f- 1 (2f(y))) 
and is continuous. Similarly we establish that ui2 '■ G — >G is a homeomorphism 
onto its image. Since / and f~ l are analytic, both w n and w^ 1 are analytic, and 
hence W is analytic. □ 

Theorem 4.2. If f : I —> K is analytic, then the graph G(f) of f is the attractor 
of an analytic interpolation IFS W — {G;wi(x,y),ui2{x,y)}, where G is a neigh- 
bourhood ofG(f). 

Proof. Let / : / — > K be analytic and have graph G(/). We will show, by explicit 
construction, that there exists an affine map T : K 2 — > K 2 of the form 

T(x, y) — (ax + h, cx + dy), ad ^ 0, 

such that T(G(/)) = G(g) is the graph of a function g : [0, 1] — > K that obeys 
the conditions of Lemma 4.1 (wherein, of course, you have to replace / by g). 
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Specifically, choose L(x) = ax+h so that T(I) = [0, 1] and then choose the constants 
c and d so that 

(x-h\ c(x-h) 
\ a J a 

satisfies the conditions (i) and (ii) in Lemma |4~T| To show that this can always be 
done suppose, without loss of generality, that I = [0, 1] so that a — 1 and h = 0, 
and let d = 1. Then to satisfy condition (i) and (ii) requires that 



max 

xG [0,1] 



/' (x/2) + c 



, max 

££[0,1] 



/' {x/2 + 1/2) 



/'(*) 



<2, 



which is true when we choose c to be sufficiently large. Finally, let W = {<G; w\, W2} 
be the IFS, provided by Lemma |4~lj whose attractor is G{g). Then W = {G;w\ = 
T~ x o wx o T, w 2 = T^ 1 o w 2 o T}, where G = T _1 (G), is an analytic interpolation 
IFS whose attractor is the graph of g. 

□ 

The following example illustrates Theorem |4.2| 

Example 4.3. (the exponential function) Consider the analytic function f(x) = 
e x restricted to the domain [1, 2]. An IFS obtained by following the proof of The- 
orem 14.21 is 



W = {R 2 ;w 1 (x,y) = {x/2 + 1/2, Jejj), w 2 (x, y) = {x/2 + l,e^y)} . 

Therefore the attractor of the analytic IFS W is the graph of / : [1,2] — > M, 
f{x) = e x . The change of coordinates (x,y) — > (y,x) yields an analytic IFS whose 
attractor is an arc of the graph of ln(x). 

5. Fractal Continuation 

This section describes a method for extending a fractal function beyond its orig- 
inal domain of definiton. Theorem 15.41 is the main result. 



Definition 5.1. If I C J are intervals on the real line and f : I —> R M 1 and 
g : J — » R M_1 , then g is called a continuation of / if / and g agree on /. 

The following notion is useful for stating the results in this section. Let X°° 
denote the set of all strings a = (Ti<J2<Jz • • • , where Cfc G T for all k. The notation 
<7i(72 • • • &m stands for the periodic string a±a 2 • ■ ■ cr m u\ ■ ■ ■ a m <j\ ■ ■ ■ . For example, 
12 = 12121 

Given an IFS W ={R M ; w n , n € 1}, if a = <JiO'2 '3 • • • G 2^°° and k is a positive 
integer, then define W a \f. by 

w a \k{x) = w ai o uv 2 o ■ ■ ■ o w CTfc (x) = {w„ 2 (. . . (w CTfc (x)) ■ • • )) 

for all x £ K M . Moreover, if each w n is invertible, define 

w 9\l ■= w e 1 low 02°- ow e^ 

Note that, in general, uw 7^ (we|i,) _1 . 

A particular type of continuation, called a fractal continuation, is defined as 
follows. Let I be an interval on the real line and let / : I —> M M_1 be a fractal 
function as described in Definition 13.41 In this section it is assumed that the IFS 
whose attractor is G{f) is invertible. Denote the inverse of w n (x,y) by 
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where F* : R M — > R M 1 is, for each n, the unique solution to 

(10) F n {L- 1 (x),F* n (x,y)) = y. 
Let 9 G 2°° and G = G(f). Define 

(11) G e]k :=w-{ k {G). 
Proposition 5.2. VKii/i notation as above, 

G c G§|i c Gg|2 c • • ■ . 
Moreover, Ggi^ is the graph of a continuous function fgi^ whose domain is 
l e \ k := L gl l(I) = L-l o L^ 1 o ... o ^(7). 

Proof. The inclusion Gg|fc_i C Ggifc is equivalent to wg k (G) C G, which follows 
from the fact that G is the attractor of W. The second statement follows from the 



form of the inverse as given in equation (10 ). □ 



It follows from Proposition 5.2 that 

oo 

(12) G e := |J Gg {k 

fe=0 

is the graph of a well-defined continuous function fg whose domain is 



h = J < 



9\k 



k=0 



Note that if 9 G X°°, then fg(x) = fg\k{x) for a; G and for any positive integer 
k. 

Definition 5.3. The function fg will be referred to as the fractal continuation 
of / with respect to 9, and {f g : 9 G 1°°} will be referred to as the set of fractal 
continuations of the fractal function /. 



Theorem 15.41 below states basic facts about the set of fractal continuations. Ac- 
cording to statement (3), the fractal continuation of an analytic function is unique, 
not depending on the string 9 G I°° . Statement (4) is of practical value, as it implies 
that stable methods, such as the chaos game algorithm, for computing numerical 
approximations to attractors, may be used to compute fractal continuations. The 
figures at the end of this section are computed in this way. 

Theorem 5.4. Let W={M. M ;w n ,n G 1} be an invertible interpolation IFS and 
let G(/), the graph of f : I —> M. M_1 , be the attractor ofYV as assured by Theorem 
\3.S\ If 9 G I°° , then the following statements hold: 

(1) 

fR if 6el°°\{l,N}, 

h = < [0:0,00) if 9 = 1, 
l(— 00, xn] if 9 = N. 

(2) fg(x) = f{x) for allxe L 

(3) IfW is an analytic IFS and f is an analytic function on I, then fg(x) — 
f{x) for all x G Ig, where f : K — > is the (unique) real analytic 
continuation of f. 
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(4) For all keN, the IFS 

has attractor Gg\}~ = G(fg\k). 



^■Wan, o w n o (w g \l)~ 



Proof. (1): Each of the affine functions L n can be determined explicitely, and it is 
easy to verify that L" 1 is an expansion. Moreover, the fixed point of L n (and hence 
also of L~ l ) lies properly between xq and Xn for all n except n = 1 and n = N. 
The fixed point of L\ is xq and the fixed point of Ln is a; at- Statement (1) now 
follows from the second part of Proposition |5.2| 

(2) : It follows from Proposition [^2] that fg{x) = f{x) for all x G I, 6 £ 1°° . 

(3) : Since F*(x,y) and L~ x are analyic for all n, each w^ 1 is analytic. Therefore 
/e(x) is analytic and agrees with f(x) on I. Hence fg(x) — f(x), the unique analytic 
continuation for x E Ig. 

(4) : It is easy to check that condition (4 I in the definition of attractor in 
Section H holds. □ 

Corollary 1. Let W ={M M ; w n , be an interpolaltion IFS, an d let G(f), the 

graph of f : I —¥ M. 1 , be the attractor ofYV as assured by Theorem 3.2 Let f be 
analytic on I and 8 G T°° . 



(1) IfW is an affine IFS, then fg(x) = f(x), where f 
analytic continuation of f . 



dM-1 



is the real 



(2) If M = 2 and W is the IFS constructed in Theorem ^.2, then fg(x) = f{x). 

Proof. For both statements, W is analytic and hence the hypotheses of statement 
(3) of Theorem 5.4 are satisfied. □ 



Remark 5.5. This is a continuation of Remark 3.3 and concerns a generalization 
of Theorem |5.4| to the case of an IFS in which the domain X of each w n : X — > X 
is a complete subspace of K M . The attractor G may not lie in the range of wg\k 
and the set-valued inverse w^jj. may map some points and sets to the empty set. 
Nonetheless, it is readily established that 

G C Gg\\ C Gg\2 C ■ • • 

is an increasing sequence of compact sets contained in . 



We can therefore define : 



G S\k = w 0\l( G ) and G e = U Gg\ k C 



k=0 



exactly as in Equations (11 



whose graph is Gg. Theorem 



and ( 12 I and the continuation fg as the function 



5.4|then holds in this setting. 



Example 5.6. This example is related to Examples 3.5 and 3.6 Let G p be the 
attractor of the affine IFS 



(]R 2 , wi(x, y) = (0.5a;, 0.5a; 4- py), 102(0;, y) — (0.5x + 1, — 0.5a; +py + 1), 



where p € (— 1, 0) U (0, 1) is a parameter. 

When p = 0.25 the attractor Go. 25 is the graph of the analytic function / 
[0, 2] -> K, where 

f(x) = x(2 - x) 
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Figure 2. See Example 5.6 




Figure 3. Added detail for part of Figure [2] showing additional 
continuations near the ends of the original function. 



and, according to statement (3) of Theorem |5 ,4| the unique continuation is fg(x) = 
x(2~x) with domain [0, oo) if 9 — 1, domain (— oo, 2] if 9 = 2, and domain (— oo, oo) 
otherwise. 

Whenp = 0.3 the attractor is the graph of a non-differentiable function and there 
are non-denumerably many distinct continuations fg : (— oo, oo) — > R. Figure 
[2] shows some of these continuations, restricted to the domain [—20,20]. More 
precisely, Figure [2| shows the graphs of fe\±{x) for all 9 € {1, 2}°°. The continuation 
fj(x), on the right in black, coincides exactly, for x G [2,4], with all continuations 
of the form fi p {x) with p e {1, 2}°°. To the right of center: the blue curve is G2111, 
the green curve is G2211, and the red curve is 62221- On the left: the lowest curve 
(part red, part blue) is G2222, the green curve is G1222, the blue curve is G1122, and 
the black curve is G1112. Also see Figure [3] 

For p — 0.8 the attractor G .s is the graph of a fractal function / .s whose graph 
has Minkowski dimension (2 — ln(5/4)/ln2). This graph G .8 is illustrated in the 
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Figure 4. See Example |5.6| The colors help to distinguish the 
different continuations of the fractal function whose graph is illus- 
trated in black near the center of the image. 

middle of Figure [f] The window for Figure [3 is [-10, 10] x [-10, 10] C M 2 , and f 
is the (unique) black object whose domain is [0, 2]. Figure [4] shows all continuations 
/ 9 |4(x)for0G{l,2}°°. 

Example 5.7. (continuation of a nowhere differentiable function of Weier- 



strasse) We continue Example 3.9 It is readily calculated that, for all £ € [0, 1), 

^l 1 ^)^) = (2a:, (y - sin27re) /£), w-^x, y) = (2x - 1, (y - sin27ra;)/£) 
from which it follows that 



k=Q 



7T.T 



with domain [0, oo) if 8 = 1, domain (— oo,2] if 6 = 2, and domain (—00,00) 
otherwise. In this example, all continuations agree, where they are defined, both 
with each other and with unique function defined by periodic extension of equation 



When W in Theorem |5.4| is affine and M — 2 write, for nel. 

F n (x, y) = c n x + d n y + e n . 

We refer to the free parameter d n £ K, constrained by \d n \ < 1, as a vertical 
scaling factor. If the verticle scaling factors are fixed and we require that the 
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Figure 5. See Example |5.8| 



attractor interpolate the data {(xi,yi)}fL , then the afhne functions F n are com- 
pletely determined. 

Example 5.8. The IFS W comprises the four afhne maps w n (x,y) = 
(L n (x), F n (x, y)) that dehne the fractal interpolation function / : [0, 1] — > R speci- 
hed by the data 

{(0, 0.25), (0.25, 0), (0.5, -0.25), (0.75, 0.5), (1, 0.25)}, 

with vertical scaling factor 0.25 on all four maps. Figure [5] illustrates the attractor, 
the graph of /, together with graphs of all continuations fijki where k, I € 
{1,2,3,4}. The window is [-10, 10] 2 . 

Example 5.9. The IFS W comprises the four afhne maps with respective vertical 
scaling factors (0.55,0.45,0.45,0.45) such that the attractor interpolates the data 

{(0,0.25), (0.25,0), (0.5,0.15), (0.75,0.6), (1,0.25)}. 

Figure [6] illustrates the attractor, an afhne fractal interpolation function / : [0, 1] — > 
K, together with all continuations fijki where i, j, k,l € {1,2, 3, 4}. The window is 
[-20, 20] 2 . 
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Figure 6. See Example [ST9 
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Figure 7. See Example 5.10 



Example 5.10. Fig ure [7| shows continuations of the fractal function / : [0, 2] — > R 
described in Example 1 3 . 7] ~These are the continuations fijki(x) (ijkl € {0, l} 4 ). The 
a;-axis between x = —10 and x — 11 and y-axis between y — —100 and y = 2 are 
also shown. The graph of / is the part of the image above the x-axis. 
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In order to describes some relationships between the continuations {fg : 9 € I 00 } 
(see the previous examples), note that, for any finite string a and any 9, 9' G I°° , 

for all x € I a . Consider the example J = [0, 1], N — 2, and 

Li (x) = - x, Lo(x) = — x H — . 

u ; 2 ' v 7 2 2 

It is easy to determine I a for various finite strings a, some of these intervals illus- 
trated in Figure [8] For example, we must have 

for all x £ I 2 2 = [—3, 1] and for all 9, 9' E but, as confirmed by examples, it 
can occur that f2i2(%) 7^ f22i(x) for some x £ I212 H I221 = [—3, 3]. 

There is a natural probability measure on the collection of continuations on 1°° 
defined by setting Pr(#i = 1) = 0.5 for all i — 1,2, .., independently. Then, because 
many continuations coincide over a given interval, we can estimate probabilities for 
the values of the continuations. For example, if N = 2, I = [0, 1], and ai — 0,2 — 
1/2, then 

Pv(f e (x) = f T (x)\xe[l,2])>l/2; 
Pi(fe(x) = Mx)\xe[l,2])>l/4; 
Pr(fe(x) = f^(x)\xe[l,2])>l/8; 

Pr(/fl(aO = /2^22l^) ' x e I 1 ' 2 D - V 2 " for all n = 1, 2, • • ■ 

n 

In this sense, Figures [2j [5j and [7] illustrate probable continuations. 
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Figure 8. An example of domains of agreements between ss ana- 
lytic continuations. See the end of Section [5] 
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6. Uniqueness of fractal continuations 

This section contains some results concerning the uniqueness of the set of fractal 
continuations. Our conjecture is that an analytic fractal function / has a unique 
set of continuations, indpependent of the particular IFS that generates the graph 
of / as the attractor. More precisely, suppose that G(f), the graph of a continuous 
function / : I — > R, is the attractor of an analytic interpolation IFS W with set of 
continuations {fg : 8 g as defined in Section|5j and the same G(/) is also the 
attractor of another analytic interpolation IFS W = {M. M ;w n ,n £ 1} with set of 
continuations {fg : 6 € 1°°}. The conjecture is that the two sets of continuations 
are equal (although they may be indexed differently). This is clearly true if / is 
itself analytic, since an analytic function has a unique analytic continuation. In this 
section we prove that the conjecture is true under certain fairly general conditions 
when / is not analytic. Recall that the relevant IFSs are of the form 

W = {X c K 2 ; w n (x,y) = (L n (x),F n (x,y)), nel}, 
L n (x) = a n x + b n 

for n = 1, 2, . . . , N. The first result concerns the extensions fj(x) and fjf(x). This 
is a special case, but introduces some key ideas. 

Theorem 6.1. Let W and W be analytic interpolation IFSs, each with the same 
attractor G(f) = G(f) but with possibly different numbers, say N and N of maps. 
Then 

/iO) = fi(x) and f w (x) = f~(x) 
for all x € K such that (x, y) 6 X for some y£R. 

Proof. As previously mentioned, it is sufficient to prove the theorem when / is not 
an analytic function. In this case / must not possess a derivative of some order at 
some point. By the self-similarity property p~| mentioned in the Introduction, / 
must possess a dense set of such points. Hence, as a consequence of the Weierstrass 
preparation theorem |11) . if a real analytic function g(x,y) vanishes on G(/), then 
g(x, y) must be identically zero. Now, since L\oh\ = L\ o L\ it follows, again from 
the self-similarity property, that [w\ owi) (x,y) = (w{ ouii) (x,y) for all (x, y) £ 
G(f). Then 

g(x, y) = (wi owx-wxo wi)(x, y) 
vanishes on G(f). Hence wi o uT\ = w\ o w\ for all {x, y) in X. It follows, on 
multiplying on the left by u>J~ and on the right by u>^ that 

w^ 1 o wi — wi o , 

and similarly that w± o uT\ 1 = w\ 1 o w±. 

Now suppose that (x,y) E G{f T ). Then (x,y) € G(/jy fe ) = w^ 1 o ... o 1 (G(/)) 
for some k. Hence we can choose I so large that 

Wi O Wi O ■ ■ ■ Wi _ 

(x, y) € W\ o wi o • • • w\ o w 1 o ... o w 1 (G(/)) 



I times 



, Wl o W\ o ■ ■ ■ W\ 

Cw^o-.-o w^ 1 » — v ' (G(/)) 

I times 



cG(/) 
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which implies 



0, y) e 



Wi 1 o wi 1 o 



.Wl 



I times 



- (G(f)) 



when I is sufficiently large. Hence G(fj) C G(/=). The opposite inclusion is proved 
similarly, as is the result for the other endpoint. □ 

6.1. Differentiability of fractal functions. We are going to need the following 
result, which is interesting in its own right, as it provides detailed information about 
analytic fractal functions. 

Theorem 6.2. Let W be an analytic interpolation IFS of the form given in equa- 



tion (13 ) with attractor G(f). Let c G [0,oo) and d G [0, 1) be real constants such 



that | J^F n (x, y) | < c and < -j^F n (x,y) < da n for all (x,y) in some neighbor- 
hood ofG(f), for all n G X. The function f : [xq,xn] — )■ K is lipschitz with lipshitz 
constant A = a c(l — d)" 1 , where a = min{a„ : n — 1, 2, N}. That is: 

(14) 1/ («)-/(*)! < A \s-t\ for alls, te[x Q ,x N }. 

Proof. Consider the sequence of iterates k = 0, 1,2, . . . 

f k+ i(x) = (Wf k ) (x) = F n (L- 1 (x),f k (L- 1 (x))) 

for x G [x n -i,x n ],n = 1,2, ...,N, where W is as defined in the statement of 
Theorem 3.2 Without loss of generality, suppose that {fk} is contained in the 
neighborhood X of G(f) mentioned in the statement of the theorem. It will first be 
shown, by induction, that f k is lipschitz. Suppose that f k {x) is lipshitz on [a;o,2Ar] 
with constant A. Then, for all s,t G [x n —i, x n ], n — 1,2,..., N we have, by the 
self-replicating property, by the mean value theorem for some C) <= ^> and by 
the induction hypothesis that 

\f k+ i(s)-f k+ i(t)\ = \F n (L-\s)J k (L- 1 (s)))-F^ 

< | F n (L- 1 (s) , f k (L- 1 (t ) ) ) - F n (L- 1 (t ) , f k (L- 1 (t ) ) ) | 
+ | F n (L~ 1 (a) , f k (L~ 1 (a) ) ) - F n (L~ 1 (a) , f k (L~ 1 (t ) ) ) | 
9 



9x 



F n {x,y) 







Fn(x,y) 



+ \f k {L-\s))-f k {L-\t)))\ 

< (s^'c + Ao; 1 d) |s-t| < (a _1 c+Ad) |s-t| 
= (a _1 c + a- 1 c(l - d) _1 d) |a-i| = A |s-i|. 

Now suppose that s < i and s G [x m _i,x m ] and t G [x„_i,a; n ] where m < n. Then 

\fk+i(s) ~ fk+i(t)\ < |/fc+i(s) - f k +i(x m )\ + \f k+ i(x m ) - f k+ i(x m+ i)\ 

+ --- + \fk+l{Xn)- fk+l{t)\ 

< A |s - x m \ + A |i m - x m +i\ H h A |x„ - £| = A \s - t\ . 



Therefore f k (x) is lipshitz on [xo,xtv] with constant A for all k. 
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Now we use the fact that {f k } converges uniformly to / on [xo,x/y], specifically 
max{|/(x) - f k (x)\ : x € [x ,x N ]} = max{| W k f(x) - W k f (x)\ : x G [x„,x N ]} 

< s k max{|/(x) - fo(x)\ : x G [x ,x N ]} 



for s as in the proof of Theorem 3.2 The unform limit of a sequence of functions 
with lipschitz constant A is a lipshitz function with constant A. □ 



For an interpolation IFS W of the form given in equation ( 13 ) with attractor 
G(f), consider the IFS 

£: = {/= [x ,x N ]; L n (x) = a n x + b n ,n G I}, 

and let 

oo 

£>£ = (J C k ({x ,xi, . . . , x N })\{x ,x N }, 

V w = {(x,f(x)) : xeV c }. 

The set 2?yy will be referred to as the set of double points of G(f). The standard 
method for addressing the points of the attractor of a contractive IFS [5] can be 
applied to draw the following conclusions. If (x, y) is a point of G(f) that is not a 
double point, then there is a unique a G T°° such that 

x = lim L a \ k (s), 

(x,y) = lim w a \ k (s,t), 

k— f oo 

where the limit is independent of (s, £) G X. If (x, y) is a double point, then there 
exist two distinct strings a such that the above equations hold. In any case, we use 
the notation 7r : I°° -> K, 

tt(ct) := lim L a \ k (s), 

k— ¥ oo 

which is independent of s G M. 

Theorem 6.3. Let W 6e an analytic interpolation IFS of the form given in equa- 
tion (13 ) with attractor G(f) and such that < \dF n (x,y)/dy\ < a n for all 
(x,y) G X and for all o£l. If x is not a double point of G(f), then f is differen- 
tiable at x. 

Proof. For now, fix k G N. Define the function Q : R fe+1 ->• R by 

Q(x!,x 2 , ...,x k ,y) = F <Tl (xi,F <T3 (x 2 , ■ ■ ■ F ak (x k ,y) ...)). 

where the domain of Q is all (iEi,X2, . . . ,x k ,y) such that (xj,y) € X for all j = 
1,2,..., k. Let 

d d 
H n (x,y) := —F n (x,y), K n (x,y) := —F n (x,y) 

for 7i G X. Because F n [x,y) is analytic for each n, there are constants c G [0,oo) 
and d„ G [0, a n ) such that 

(15) \H n {x,y)\ < c, \K n (x,y)\ < d n 
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for all n £l and for all (x,y) in some neighborhood of G(f). Using the notation 



Q xl (xi,X2, ...,x k ,y) ■ 



Q y (x 1 ,x 2 , ■ ■ ■ ,x k ,y) ■ 



_d_ 
dxi 
d_ 
dy 



Q(xi,x 2 , ■ ■ -,Xk,y), 



Q(xi,x 2 , ...,x k ,y) 



we have the following for I = 1, 2, . . . , k: 

Q Xl (x 1 ,x 2 , ...,x k ,y) = K ai (xi, F CT2 (x 2) . . . F ak (x k ,y) ■ ■■))■ K a2 (x 2 , F<r 3 (x 3 , . . . F ah (x k) y) ...))■... 

K ai _ 1 (xi- 1 , . ..F ak (x k ,y) ...))• H ai (xi,F ai+x (xi + i, . ..F ak (x k ,y) ■■■)); 



Q y (xi,x 2 , ...,x k ,y) = K (T1 (x 1 ,F (T2 (x 2 ,. . .F ak (x k ,y) ...))■ K a2 (x 2 , F a3 (x 3 , . . . F ah {x k ,y) ...))■... 
K <Jk _ l (x k - 1 ,F <Jk (x k ,y)) ■ K„ k (x k ,y). 

Using the intermediate value theorem repeatedly we have, for some rj l £ [xi,xi+8xi], 
I = 1,2,..., k, and £ £ [y,y + Sy]: 

AQ :=Q(x 1 + 8xx,x 2 + 8x 2 , ...,x k + 5x k ,y + Sy) - Q(x 1 ,x 2 , ...,x k ,y) 

= [Q(xi + Sxi, x 2 + 6x 2 , . . . ,Xk + Sx k ,y + Sy) - Q(x 1 ,x 2 + Sx 2 , . . . ,x k + 5x k , y + Sy)]+ 

[Q{xi,x 2 + fa 2 , ...,Xk + Sx k ,y + Sy) - Q(x 1} x 2 ,x 3 + Sx 3 , ...,x k + Sx k ,y + Sy)} + ... 

[Q(xx,x 2 , . ..Xk-i,x k + Sx kl y + Sy) - Q{x 1} x 2 , ...,x k ,y + Sy)} + 

[Q(x 1 ,x 2 , ...,x k ,y + Sy)- Q(x x ,x 2 , . . .,x k ,y)} 
= Q Xl (rj 1 ,x 2 + Sx 2 , ...,x k + Sx k ,y + Sy)Sxi + 

Q X2 (xi,r/ 2 ,x 3 + Sx 3 ...,x k + Sx k ,y + Sy)Sx 2 H h 

Qx k (xi,x 2 , x 3 ... x k -i,r] k ,y + Sy)Sx k + 

Q y (x 1 ,x 2 ,x 3 . . .,x k ,£,)Sy. 

Let a — tt^ 1 (x), which is well defined since x is not a double point, and 




fk(x) = F ai (x 1 ,F cr2 (x 2 , ...F ah (xk,fo(x k )) ■■■)). 



Fix k and a and let both x and (x + Sx) lie in L a \ k ([0, 1]). Define 



■j = (io-y) 1 (x) 
y = .f((L alk )-\x)) 



Xj + Sxj = (L a \j) 1 {x + Sx) = (L a y) 1 (x) + (a, 
y + Sy = f((L alk )- 1 (x + Sx)). 



...a aj ) l 6x 
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Then 

[f(x + 6x)-f(x)]/6x 

= [F ai (xi + 5x 1 ,F lT2 (x 2 + 5x 2 , ...F ak (x k + 5x k , y + Sy)... ))- 

Fa i (xi,F a2 (x 2 , ■ ■ ■ F a k (x k ,y)... ))]/5x 
= AQ/Sx 

= [QxAVi,X2 + 6x 2 ,.. .,x k + Sx k ,y + 5y) ■ (a CTl ) _1 + 

Qx 2 (xi,V2i x 3 + 5x 3 . . . ,x k + 5x k ,y + Sy) ■ (a^a^) -1 + . . . 
Qx k (xi,x 2 ,x 3 . . .x k -i,r] k ,y + 6y) ■ ( a ai . . . a ak 
Q y (x 1 ,x 2 ,x 3 ...,Xk,Q-{a 

) _1 • 5x) — f((L a \ k )~ (x)) /((a CTi • • • a ak ) ■ 5x) 

and 

[f(x + 6x)-f(x)]/5x 

= [cia~l + c 2 d 2 ^a~\a~\ + c^^.^la^a^l + ... 

+ c k d k _xd k ^ 2 . . . rffe ; (fe_i)(a (Tl . . . a ak ) 1 ] 
- _£L -i- ^1. ( hiL c 2 d 3 ^d 3 . 2 ^ c k d k ,xd k _ 2 . . . <ifc,(fc-i) 

/((^ |fc )- 1 ( 2; + <5x))-/((^ |fc )- 1 (x)) 



where 



Qj/(^l,^2,^3 ■ ■ ■ ,Xk,Q 

(Iq- 1 d(j2 ^(73 • • * ^CTfc 



ci = H r7l (r] 1 ,F <J2 (x 2 + 5x 2 , ■ ■ ■ F ak (x k + 5x k ,y + Sy) . . . )), 
c 2 = FLfj 2 {r]2,F a3 (x 3 + 5x 3 , . . . F ak (x k + 5x k , y + 6y) . . . )), 

Ck = H ak (i] k ,y + 5y), 

and 

cfe,i = K ai (xi , F a2 (r? 2 , F CT3 (ar 3 + &r 3 , . . . (x k + 5x k ,y) ■■■))) 

d 3 ,i = K ai (xi, F a2 (x 2 ,F a3 (rj 3 ,F a3 (x 4 + 5x 4 , ■ ..F t7k {x k + 6x k , y) ■ ■ ■ )))) 

d.3,2 = K 02 [x 2 , Ffj 3 (r? 3 , F a4 (x 4 + 5x 4 , . . . F ak (x k + Sx k ,y) . . .))) 

di,i = K ai {x\ , F a2 (x 2 ,... F ai _ 1 [xi—\, F ( j l (r] l , F a[+1 (x i+1 + 6x l+1 , ...F ak (x k + 5x k ,y) 
di2 = K a2 {x 2 , F a2 (x 2 ,...F ai _ 1 {xi—\, F ( j l F ai+1 (x i+1 + 5x l+1 , . . . F ak (x k + 5x k ,y) 

di,i-i = F ai _ 1 {xi—\, F (7l (rii,F ai+1 (x l+ i + 5xi +1 , ...F ak (x k + 5x k , y) . . . ))) 
4,i = K ai (x u F„ 2 (x 2 ,... F ai _ 1 {xi-\,F„ t {xi,F ai+1 (x l+1 , . . . F CTfc (r? fe , y) . . . )))) 



d k ,k-i = K ak _ 1 (x k ^ 1 ,F ak (i lk ,y)). 
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Note that the Ci and di ym depend explicitly on k (which so far is fixed). It follows 
that 



f(x + Sx)-f(x) ( Cl | d 2 ,i C 2 | ^3,1^3, 2 c 3 



5x 



dksdka ■ ■ ■ dk,(k-i) Ck 



L (7\ Uj (72 ' ' ' "'Cfc — 1 "'Cfc 



< 



Q y {x 1 ,x 2 ,x 3 ... ,x k ,£) 



< X 



y (a;i,a;2,X3 . . . ,x fe ,C) 



/(S fc (a: + fa))-/(5 fc (j)) 
S*(z + fe)-Sfc(a:) 



^(Ti Ci(T2 ^*C"3 * * ' ^CTfc 



the last inequality by Theorem 6.2 The above is true for all x, (x + dx) € [xq, Xff], 
5x ^ 0, k e N. We also have 



(^1,^2,^3 ■ ■ ■ ,X kl Q 



\K ai {x 1 ,F c , 2 (x 2 ,...F ah (x k ,£)...))\ 



I K a2 (x 2 , F„ 3 (x 3 , . . . F ak (x k , • • • ) ) I 



<n^<^ 



for some C € [0, 1), the last inequality by equation (15 ). Hence, for any e > 0, we 
can choose k so large that 



(16) 



f(x + Sx) - f(x) ^ c m "t-t d, 



5x 



k,l 



m—1 



1 = 1 



< e/3. 



Note that, by their definitions, for fixed x, the c m s and d kf is depend upon both k 
and Sx. Our next goal is to remove the dependence on both k and Sx. For all I and 
all k > I define 

(17) C ai :=H ai (xi,f(xi)=H ai (xi,F ai+1 (xi + i,...F ak (xk,f(xk))---)) 
D ai := K ai (xi,f(xi) = K rJl (xi,F ai+1 (xi +1 , . . . F ak (x k , f(x k )) ...)). 

We are going to show that, for all e > and for Sx sufficiently small, 



(18) 

and that 
(19) 



' C. Tr 1 D 



1 a " m , , a, 

m— 1 7,1 Z=l 



A: m—1 , 

Cm TT «/c,i 



1:111: ii 



i a ° m 7 , a o 

m—1 m l—l 



oo ^ m—1 j-. fe ^ m-1 n 



m—1 m / — 1 



Z=l 



<e/3, 



< e/3, 



which taken together with inequality ( 16 ) imply 
(20) 



f(x + Sx)-f(x) 



5.x 



m=l °™ 1 = 1 



< e. 
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For \8x\ sufficiently small 



k , i 



a„ 



_ ~ n 



< 



=1 

k 

E 

m— 1 



1 = 1 



'k.l 



\c„ 



n 



I A, 



E^n 



ii>„ 



l k,l\ 



< e/3. 



z=i 



The last inequality above follow from, for fixed fc, the continuous dependence of 
the c m s and dk ;S on their independent variables, and comparing C CTwi with c m 
and Dcr, with ci^; using the equalities ( |17 [ )■ (We need \6x\ small enough that 



+ Sx lies in L CT i fc ([0, 1]).) We have established (18 I. Concerning inequality 



(19 ), by equation (15 ) the c m 's are uniformly bounded and, for some (x,y), 



we have \d k ,i\ = \K ai (x,y)\ < d ai < a ai . Therefore \d k ,i/a ai \ < \d a Ja ai \ < K for 



some constant K < 1. So inequality (|19 J foll ows from the absolute convergence of 

oo — m— 1 



the series 



n 

1=1 



D„ 



From Equation (20 



it follows that 



m = E 



c_ 



n 

n — 



Z = l 



□ 



Note that the last equality in the above proof actually provides a formula for the 
derivative at each point that is not a double point. 

6.2. Unicity Theorem. We conjecture that the uniqueness of the set of contin- 
uations holds in general. The following theorem provides a proof in M 2 under the 
assumption that the derivative f'(x) does not exist at all points x, although we 
conjecture that uniqueness holds in M. M , M > 2, and it is sufficient to assume that 
f(x) is not analytic. It is also assumed that there is a bound \dF n (x,y)/dy\ < a n , 
where the a r , 



are as given in equation ( 13 



As an example, consider the case of 
affne fractal interpolation functions, where F n (x, y) — (a n x + b n , c n x + d n y + g n ). 



Then for Theorem 6.4 to apply we need \d n \ < a n for all n. 

{X C Rj, w n {x, y) = (L n (x),F n (x,y)), n G 2} and 
(L n (x),F n (x,y)), w£l} be analytic interpolation IFSs 

< a n 



Theorem 6.4. Let W 

W = {XcR 2 ;w IL (x,y)) 
as in equation 



13 ) such thatO < \dF n (x,y)/dy\ < a n andO < dF n (x,y)/dy 

for all (x, y) £ X, for all n £ X. If both W and W have the same attractor G(f) 
such that f'(x) does not exist at x — x n , for all n = 0, 1, 2, ... , N, then VV =W. 

Proof. For simplicity we restrict the proof to the case N = 2. The proof of the 
result for arbitrary many interpolation points is similar. 

We first prove that the set of double points of G(f) with repect to W is the same 
as the set of double points of W. The interpolation points for W are {0, x\ , 1} and 
the interpolation points for W are {0, x\, 1}. By Theorem 1 



6.3 



f{x) is differentiable 

at all points that are not double points with respect to W and also at all points 
that are not double points with respect to W. Moreover, f{x) is not differentiable 
at all double points with respect to W and also not differentiable at all points 
which are double points with respect to VV. (Otherwise f{x) must be differentiable 



FRACTAL CONTINUATION 



2:S 



at x\ which would imply that f(x) is differentiable everywhere, contrary to the 
assumptions of the theorem.) It follows that f(x) is not differentiable at x if and 
only if a; is a double points with respect to W if and only if x is a double point with 
respect to W. 

We next prove that w n (x, y) = w n (x, y) for all (x, y) G G(f) and n = 1, 2, . . . , N. 
Since X\ is a double point of G(f) with respect to W there must be a\k ^ such 
that w a if.(xi, f{xi)) = (xi, f(xi)). Since Xi is a double point of G(f) with respect 

to W there must be a\k such that w~^(xi, f(xi)) = (xi, /(xi)). It follows that 
™(Z\k)( w (<r\k)(xi, f(xi))) = (xi,f(xi). Since W(g\k) ° w (v\k) ■ G(f) G(/), we can 
write o w^ik) (x, y) — w(x,y) = (L(x), F(x,y)) where, similar in form to the 

functions y) and w n {x, y) that comprise the two IFSs, L{x) = ax + h is a real 

afhnc contraction and F(x,y) is analytic in a neighborhood of G(/) and has the 

property, by the chain rule, that ^(x,y)} < a in a neighborhood of G(f). It is 

also the case that axi + h = x\ and F(L (x), f(L X (x))) = f(x) in a neighborhood 
of xi and £(xi) = x\, F{x\, f(xi)) = f{x\). Using the analyticity of F(x,y) in x 
and y, 

f( Xl +6x)~ f( Xl ) = F{L- 1 {x 1 + Sx), f{L-\x 1 + 5x))) - T(L^ 1 (xi), f(L^ 1 (xi))) 
5x Sx 
= F x (x!, fixi))^ 1 

.-nt t, nn— l - f(xi)) fx s 
+ ™y(zi,/0Ei))a + o(ojc). 

This implies that the following limit exists: 

( /(x x + Sx) - /On) 1 (/(a; 1 +Q- 1 fa)-/( a ; 1 )) | 

lim < Fy(xi,f(xi))a } 

Sx^o 1 Sx a L Sx j 

= (1 - F y (x u fixi^cT^fixt) = F^xufixx))^ 1 , 

which implies 

(a- F 1/ (a;i,/(ar 1 ))) 

We have shown that if <r|fc 7^ then /(x) is differentiable at x\, which is not 
true. Therefore o\k — which implies x\ — x\ and hence w n (x,y) — w n (x,y) for 
a dense set of points (x, y) on G{f). It follows that w n {x,y) — w n {x,y) for all 
(x,y) G G(/) and n = l,2. 

To show that W =W, i.e., that w n (x,y) — w n (x,y) for all (x,y) G X, define an 
analytic function of two variables, a : X — > R by a(x,y) := w n (x,y) — w n (x,y) for 
all (x, y) E X. It was shown above that a(x,y) — for all (x, y) G G(f). That 
a(x,y) = for all (x, y) G X follows from the Weierstass preparation theorem 

ca. □ 
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